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Abstract 

We compute generating functions for number of U(N)(SU(N)) singlets in Fock space 



in several space dimensions. The motivation to find the explicit form of the functions is 
from the numerical approach to supersymmetric Yang-Mills quantum mechanics, based 
on Fock space. Incidently the functions give many important insights into the quantum 
mechanical models based U(N)(SU(N)) gauge group. 

On 

1 Introduction 

o 

Yang Mills theories in zero volume limit provide a natural arena where quantum mechanics 
with the singlet constraint emerge due to the Gauss law The gauge fields Af, where i 
and a are spatial and color indices respectively, become now the coordinate operators xf 
in emerging Yang-Mills quantum mechanics (YMQM). The resulting system is extremely 
difficult to solve and accordingly there are no exact solutions in the literature. 

The model and its supersymmetric extension (SYMQM) play an important role in 
quantum mechanical description of the membrane [2] and the supermembrane [3], i.e. 
they give a regularized description of a (super)membrane. The D = 9 + 1 dimensional 
SYMQM model [3] became famous due to the BFSS conjecture [5], relating the later to 
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M-theory, and since then has been studied in number of papers ( we refer to the existing 
reviews [6] ). However there are no nonperturbative calculations of e.g. the spectrum of 
BFSS matrix model because of the high complexity ( of e.g. the Fock space ) that makes 
any numerical approach difficult. 

A particularly useful approach is the cutoff method (Jj, which consists of representing 
the hamiltonian in truncated Fock space and diagonalizing the resulting finite matrix. In 
the number of papers [TJ EJ EJ [101 E] the method, applied to the D = 2,4, SU(2) system, 
proved to be very fruitful giving the nonperturbative values of the spectra, wave functions, 
the Witten index etc. The Fock space approach can also be applied to the systems with 
large number of colors N |12l \13\ [Tl"l fT5j [T6] where the nonperturbative results may shed 
some new light into the area of large N quantum field theory. 

An important step in applying the cutoff method is the construction of the SU(N) in- 
variant basis in Fock space. This however proved to be very time consuming in some cases 
[17) . We focus in this paper on computing the number of SU(N) singlets in Fock space 
analytically hence facilitating the numerical considerations. The difficulty in computing 
this number arises due to many identities ( in fact an infinite number of them, emerging 
from the Cayley-Hamilton theorem ) among SU(N) singlets that make them linearly de- 
pendent. There are, of course, very well known algorithms to produce the singlet states 
[18J. One simply takes the trace operators, e.g. Tr(xiXj . . .), X{ = xfT a where T a 's are 
SU(N) generators in the fundamental representation, and act with an arbitrary number 
of products of trace operators on e.g. the Fock vacuum. Such set of states certainly 
spans the whole Fock space however they are in general linearly dependent. One can use 
the Cayley-Hamilton theorem to find the linearly independent ones but this is in general 
tedious when the number of X{ matrices becomes large. 

There is however a way to circumvent this problem when N = oo since, in that case, 
there are no additional identities among singlets. Furthermore, if one is interested only 
in single trace states then the famous Polya theorem [19] can be applied to compute 
the number of singlets explicitly. The first computation of this kind is by Sundborg [20J, 
Polyakov [21] and Aharony et al. [22] in the context of weakly coupled Yang Mills theories 
and by Semenoff et al. [23] in the context of BFSS matrix model. 
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In quantum mechanical systems the property, that the single trace states are the most 
important ones in large N limit, is not always valid. One can verify [12] that in the 
case of the anharmonic oscillator the single trace states do not reproduce the correct [21] 
large N spectrum. A different example is the free hamiltonian [25] where it was explicitly 
shown that all the 1/N k terms play an important role in the large N limit due to the 
distinguished role of the bilinear operators Tr(x 2 ). 

In this paper we neither assume that N is large nor do we consider only the single 
trace states. The method that we will exploit in details in sections 2 and 3 is the very well 
known character method in group theorjQ. In this framework we compute the generating 
functions for numbers of singlets in sectors with given number of bosonic/fermionic quanta 
and therefore give the exact number of singlets. In section 2 we discuss only the D = 2 
SYMQM. The generating function of number of gauge singlets can be related, in this case, 
to the Witten index of particular models. In this way we obtain some physical constraints, 
on seemingly unrelated group theory numbers, which we find very interesting. In section 
3 and 4 we do the analogous computation only for D = 4, 6, 10 and relate our generating 
functions with the partition function of the weakly coupled Yang- Mills theories on R x S . 

While this work was done a related paper [2H] by F. Dolan was published. There is 
a significant overlap between [2S] and our paper although the motivation as well as the 
detailed calculations are in different spirit. 



2 The character method for D = 1 + 1 SYMQM 

Let us consider the most general state in D = 1 + 1 SYMQM 

|,)=T fel .., nflCl ... c „ F at\..at b - / t Cl .../t c -| ). 

Here | 0} is the Fock vacuum while or b and f* b are bosonic and fermionic creation 
operators being in the adjoint representation of the Lie group G. The operators obey the 
commutation and anticommutation relations [a b ,a^ c ] = 5 bc , {f b ,f^ c } = S bc . The tensor 
Tb 1 ...b ng ci...c n i s a group invariant tensor so that the state | s) is a group singlet. The 
1 l thank R. Janik for the idea. 
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number of bosonic quanta in | s) is ub and the number of fermionic quanta is np. In this 
case we say that | s) is in the (ne,np) sector. 

For the Lie group G we take U(N) ( or SU(N) ). The adjoint representation of U(N) 
will be denoted by R hence the state | s) is in the representation Sym(<&™E 1 R)x Alt{®™E-±R) 
where Alt and Sym stand for the anti-symmetrization and symmetrization of the tensor 
product respectively. Let us denote the number of U(N) singlets in (nB,np) sector by 
Dng N J F . Using the orthogonality of group characters, we have 

d"S f = J d mN)X % B 2(R)x [a ^{R\ 

with the group invariant measure d/J,uw\. The symmetric and antisymmetric powers of R, 
Xsh/m(R)> X^~aii}(.R) an d the characters x can be readily constructed using the Frobenius 
formula ( the complete calculation is given in the Appendix A ). 

The direct evaluation of Dn^n F is difficult, however we can evaluate it indirectly 
introducing the following generating functions 



oo oo 

G u ^(a,b) =J2 E D ni N l* nB (-b) nF , I a |< 1, beR. (1) 

ng=0 np=0 

The sum over is in fact finite ( due to the Pauli principle, i.e. X [ Ai t (R) = for k > N 2 
for U(N) ) but it is more convenient to work with infinite sum. In Appendix A we derive 
the following integral representation 

A similar expression was derived by Skagerstam [27J in the context of a singlet ideal gas 
in flat space. 

Above integrals can be calculated explicitly for 6 = 0. We have ( see Appendix B ) 

N 1 oo 

G u ^(a,0) = ]J JY^k) = E * nB PN(n B ), (3) 

fc=l ^ ' n B =0 

where pjv(ns) is the number of partitions of ub into numbers 1, 2 . . . , N, i.e. the number 
of natural solutions of the equation Y2k=l = n B- 
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In short the calculation goes as follows. First we use the Cauchy determinant formula 
and rewrite the expression under the integral ([2]) as a determinant of a certain matrix. 
Then we express the determinant as the sum over cycles and integrate each cycle sepa- 
rately. Finally, we use the Cayley identity to write the result in a compact form. 

According to Eqn. ((2j), G u{N \a,b) is a polynomial in variable b. The coefficient next 
to 6° is equal G u ^(a,0). We now write © as 



/ N \ N 2 

<^c»>-(np^,)£< 

\jfe=l v > J i=0 



lyVc^ia), (4) 



where c^ N \a) are polynomials in variable a and Cq^ N \o) = 1. |^ 

The case of SU(N) group is analogous. We have ( see Appendix B ) 

(N \ N 2 -l 

Ujrh)) E (-')'"' " <K, ("). 
/c — 2 / i — 

N 1 oo 

k=2 K ' n B =0 



where ?jv(ib) is the number of partitions of ub into numbers 2, 3, . . . ,N and c? (a) 
are polynomials in variable a. 

The determination of (a) or (a) is difficult for arbitrary N, however we can 

can get an idea about their structure by simply counting the constructed gauge singlets. 

2.1 Explicit construction of gauge invariant states 

The formula ([3|) can be derived by directly counting the U(N) singlets. If we introduce 
the matrices = a^Tf, where TJ, are U(N) generators □ then all the singlets are linear 
combination of the following states 



h,i 2 ,...,i N ) = (a^(a^y 2 ...^ N y»\0), (5) 



2 At this stage it is not evident that one can factor out the term G U ( N ^ (a, 0) from G U ( N ^ (a, 6) leaving 

Cj (a)'s in the form of polynomials. However it is indeed the case as we shall see in the next subsection. 
3 We use the following conventions for U(N) and SU(N) generators 

T a Tb — 2JY^ ab ^ (d a bc + ifabc)T c , 

where d a b c and f a bc are corresponding structure tensors. 
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where we used the shorthand notation (A) = Tr(A) where A is an arbitrary matrix. 
Higher powers of do not appear due to the Cayley-Hamilton theorem. Vectors ([!]) are 
already linearly independent therefore they form a basis in the space of U(N) singlets. 
For a given number of quanta ub the number of vectors ([5]) is equal to the number of 
natural solutions of the equation Y^k=i^k = ^b, i>e. there are pn&b) such vectors. 
Therefore the generating function is exactly ([3]). 

The advantage of the explicit integration over characters is that we can compute 
the number of singlets in all fermionic sectors. The direct construction of such states 
is possible introducing the matrices /t = f^T^. It follows that all the singlets can be 
obtained by acting with (a) 1 f^ 3 . . .) operators on | 0). There are however new identities 
among such vectors that make many of them linearly dependent. The process of choosing 
the independent ones is in general tedious therefore the polynomials c^ N ^ (a) are difficult 
to obtain by simply counting the constructed singlets. On the contrary, using the character 
method this is straightforward ( for given N). 

From @ we see that 

G u ^(a,b) = b N2 G u W(a,l/b), 

therefore 

cf (A °(«) =cE(«), 

hence 

D U(N) _ D U(N) (6) 

Equation ([6]) is the simplest example of constraints on Dn^n F ■ It follows that the number 
of gauge invariant states in (riB,nj?) sector is equal to the number of gauge invariant states 
in (riB,N 2 — np) sector. The identity is related to the particle-hole symmetry, i.e. the 
invariance of the system under the transformation f* a — * f a , f a — > /' a (the empty 
fermionic states are replaced by filled fermionic ones). 

Examples of c^ N \a) and cf U ^ N \a) for N = 2,3,4 are in the Appendix C. We can 
determine these polynomials for low values of N by explicitly constructing the singlet 
states. Let us start with the case of SU(2). According to character integrals we have ( 
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see Appendix C ) 

SU(2) _ 1 5(7(2) _ 5(7(2) _ 517(2) _ 1 

Cg — i, — a, c 2 — a, c 3 — i, 

and the corresponding states are 

| i) = (at 2 )* | 0), (at/t) | i), (at/t/t) | j), (/t/t/t) | 

There are no other, independent singlet states which (for SU(2) group) can be seen from 
the following simple argument. The SU(2) tensor T a ___b... in ([2]) can only be made out of 
linear combinations of products and contractions between the SU(2) primite tensors e a b c 
and S a b- However products ( and in particular contractions) of e a fe c 's can be expressed as 
linear combinations of products of Sab's. Therefore, T a ___b... is either a product of 5 a b's or 
a product of 6 a b's and one e a b c - Since for SU(2) group e a b c oc (T a TbT c ) we arrive at the 
following states 

| i) = (ot 2 ) 4 | 0), (at/t) | i)> ^t/t/t) | i)> (/t/t/t) | i)s (a t/t )(a t/t/t) | ^ 

(according to the fermion exclusion principle we have (/t/t) = and (at/t) 2 = hence 
there are no operators of this form). The last state is in fact proportional to the last but 
one. To see this we write 

(at/t)(at/t/t) oc a^t a t/t/t + a t/t a t/t/t + a t/t a t/t/t K ( a t a t) (/t/t/t) 

The above arguments are difficult to generalize for SU(N > 2) due to the additional 
completely symmetric tensors d a b c that have to be considered and the fact that the prod- 
ucts of f a bc tensors cannot be expressed in terms of products of 5 a bS. It is more convenient 
use another method, which takes advantage of the Cayley-Hamilton theorem and which 
can be generalized for SU(N > 2). In section 2.1.1. we give the details of this approach 
for the SU(3) group. 

The case of U(2) is more complicated. We have 

c^ = l, cf (2) = l + a, c U 2 {2) = 2a, = 1 + a, c?< 2 > = 1, 

and the corresponding singlets are 

|i,j) = (ot)V 2 ) J '|0)«c^ (2) = l I 



(/ f )|i,i>, (at/t)|i,i)^ c ^ 2 ) = l + a, 

I ij), I *,i> — <? (2) = 2a, 

(/W) I i,i>, I i,J> — ^ (2) = l + o, 

(/ f )(/W) — cj (2) = i. 

The only difference between the SU{2) case are the additional operators (at) and (/t) 
that have to be considered when constructing the singlet state. 

The construction of independent states for SU(3) or U(3), in all fermion sectors, is 
already very nontrivial and will be discussed in section 2.1.1. However, for np = 0, 1 we 
can construct them for arbitrary U(N) and SU(N). They are 

I ■ - \ U(N) .. 

| ll,l 2 , ■■■,IN) < * C Q = 1, 

(/W*) I ii,i 2 ,...,i7v), A; = 0,...,iV-l <— ► cf (7V) = l + a + .-. + a^" 1 , 



In- • \ S!7(JV) 1 

(/V^lO.ia,...,^), fc = l,...,JV-l <— cf^ =a + ... + o JV - 1 . 

From the explicit results for N = 2,3,4 we see that there is a relation between the 
generating functions for U(N) and SU(N) namely 

G U(N) M = ^zb G SU(N) iaib) 
1 — a 

It can be understood in terms of gauge singlets just constructed. For the case of U(N) we 
have additional trace operators (a') and (f>) ( compared to the SU(N) case) therefore 
to obtain all the U(N) singlets we have to multiply the SU(N) singlets by (a^) k , k > 
and by (P) k , k = 0,1. In terms of generating functions this corresponds to multiplying 
G su W(a,b) by 

2.1.1 The construction of SU(3) invariant states 

The explicit evaluation of the integral ^§ for N = 3 gives 

c SU(3) = x c SU(3) =a + a2 ^ c SU(3 } =a + a 2 + ^ 
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cf (3) = 1 + a + 2a 2 + 3a 3 + a\ cf ® = 2a + 4a 2 + 2a 3 + 2a 4 , cf ( 3 > = c 8 f (3) . 

It follows that D nB )n F with np > can be obtained by linear combinations of D nB = 
Q3 {kb)' According to the above results we have 

=q 3 (n B -l)+q 3 (n B -2), 
D T B f = Q3(n B - 1) + qs(n B ~ 2) + 2g 3 (n B - 3), 

SU(3) 

D n B ,3 = <l3{n B ) + a 3 (n B - 1) + 2q 3 (n B - 2) + 303(715 - 3) + q 3 (n B - 4), 

D nsA ] = 2( i^ n B - 1) + 4g 3 (n B - 2) + 2a 3 (n B - 3) + 2g 3 (n B - 4). (7) 

We now construct the states in Fock space corresponding to these numbers. The con- 
struction is done separately in sectors with given number of fermionic quanta. 
The bases for tif = and np = 1 sectors are 

I ij) = (aWjVataty" | 0) < — > 1, (8) 

and 

(a f / f ) I hj) < > a, 

(oW)|i 1 j)^o J . (9) 

The number of vectors ([8]) with given number of quanta rig is exactly q3(n B ) and the 
number of vectors ([9j) with given number of quanta n B is precisely q 3 (n B — 1) + q 3 (n B — 2) , 
i.e. the number of states (a* f') \ with n B bosons is 03(715 — 1) and the number of 
states (af p p) \ i,j) with n B bosons is 03(71/3 — 2). This result agrees with D^ff and 
D S n U B f in (0). 

Our strategy to determine the basis in sectors with np > 1 is the following. First 
we list all possible trace operators with given number of fermions. Then, with use of the 
Cayley-Hamilton theorem and some symmetry arguments, we choose only the linearly 
independent ones. The resulting trace operators acting on states | are our candidates 
for the basis. In order to find out whether they really form a basis we compare the 
number of singlets constructed in this way with Dng,np in © (or equivalently associate 
the corresponding polynomial c^ U ^ 3 \a)). 



In rip = 2 sector the possible trace operators are 



(aW/t/t), (ot/tot/t) = 0, (a 



t / t )(a t / t ) = 0, n B = 2 



(otoVV/t), (at/t)( a t a t/t) 



= 3 



(aV/WaV 1 ) = 0, (aW/tXaV/t) = 0, n B = 4. 



The other trace operators involve at least one a) hence they are linearly dependent. We 
have already taken into consideration the cyclicity of the trace. Therefore, there are four 
families of SU (3) invariant states in this sector 



(/ f3 ), n B = 0, 
(at/t 3 ), n B = 1, 
(at 2 / t3 ), (aW/t 2 ), ( a t / t )(a t / t 2 )) nfi = 2i 
(at 2 / f aV t2 ), (at 2 /t 2 at/t), ( a t / t a t / t a t/t) j ^f)(a^f 2 ), ( a t 2 /t)( a t/t 2 ), n B = 3, 
(at 2 /t a t 2 /t 2 ), (at 2 /t a t/t a t / t )j ( a t/t)(at 2 /t a t/t), (atat/t)( a taW), n B = 4, 
(at 2 /t a t 2 /t a t/t), (at 2 /t)( a t 2 /t a t/t), n B = 5, 
(a+W'/WV), n B = 6. 



Again, other trace operators involve at least one hence they are linearly dependent. 
From the 17 operators listed above only 8 are linearly independent. The linear dependence 
is due to the identity for SU (3) generators Q 



(a+ot/t/t)|i,j) + 



(10) 



In the sector with = 3 the possible trace operators are 



T( a TbT c ) = b{abT c ) + 4d abc l. 



4 The complete symmetrization over indices is 




without the conventional -k factor. 
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This identity is equivalent to the Cayley-Hamilton theorem for 3x3 traceless matrices. 
If we contract its left and right hand side with a) a a\f \ we will obtain an analogue of 
Cayley-Hamilton theorem (with fermionic matrices) 

at a t/t + /t a t a t + a t/t a t = i( a t Q t)/t + ( a t/t) a t + ( a t a t/t). 

All the linearly dependent trace operators can be derived from the last identity ( see 
Appendix D ) and one finds that there are only 8 families of independent vectors. They 



arc 



(/ t3 )Ki) — i, 

f / t3 ) I hj) < ► a, 

(atat/t 3 ) | i,j), (at/t)( a t/t 2 ) | i,j) <_> 2a 2 , 
(«t 2 / f V/t) | (at 2 / f )(at/t 2 ) | ij), (at/t )(a t 2 / t 2 ) | 3^ 

( a t 2 /t)( a t 2 /t 2 ) | ^) ^ a*. (11) 

The up = A sector is the most complicated one. There are 52 trace operators that one 
can construct however only 10 of them are independent. They are 

{a)f^)\t,j), ( a t/t)(/t 3 )| iji )^2a, 
(at 2 /t 4 ) | ij), («t 2 /t)(/t 3 ) | ij), (at/t)( a t/t 3 ) | iJ)t ( a +/t 2 )(a+/t 2 ) | i, j) <— 4a 2 , 
(at 2 /t)( a t/t 3 ) | ( a t/t 2 )( a t 2 /t 2 ) | 2a 3 , 

(at/t)( a t 2 /t)( a t/t 2 ) | ( a t/t 2 )(at 2 /t a t / t) | 2a \ 

The method to extract these 10 independent ones is the same as in the uf = 3 case. 

2.2 Physical constraints on Dng$ F 

There is a class of identities analogous to © which also have a physical interpretation. 
For example if we put b = 1 in ([2]) we obtain 

G c/ W(a,l)=0, 

therefore 

V = V D U W (12) 

np even n F odd 
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The above equation is suggestive of supersymmetry since the number of fermionic states 
( the sum over odd np ) equals to the number of bosonic degrees of freedom (the sum 
over even np ). However, since we did not specify any hamiltonian, Eq. (I12p tells us that 
the bosonic and fermionic states match, as they should, i.e. the Hilbert space of singlets 
is already "prepared" for supersymmetry. 

An interesting observation is given by equation (pQ), when a = b. The result is of a 
form of the Witten index. Indeed, if we consider supersymmetric harmonic oscillator 

H = {Q, Q 1 "} = a) b a b + p b f b , Q = a^ b f b , 

then the energy is proportional to the number of quanta hence there is only one vacuum 
state and naturally it is the Fock vacuum. Therefore, the Witten index is 1. We con- 
firm that by explicitly computing G u ^ N \a,a) using Eqn. ([2]) from which the identity 
G u( - N \a,a) = 1 follows. We also note that if we put b = a in ([T]) we obtain 

riB, np k>0 n B +n F =k 

therefore 

E (-ir «l = o- 

nB+nF=const.>Q 

The origin of the above identities lies in the dynamics of the supersymmetric harmonic 
oscillator although it is perhaps not evident at first sight. 

Another example of a hamiltonian which brings physical meaning to some identities 
including Du^n^s is the hamiltonian given by the supercharge [12] 

Q = f b cJ h + gd ahc a) V b f c . 

One can show [TS] that the in the limit of strong 't Hooft coupling A = Ng 2 — > oo the 
energies are proportional to + 2np and that the supercharges act in the subspace 
of vectors such that n# + 2np is fixed. Therefore the contribution to the Witten index 
in terms of generating function is now This quantity is not necessarily a 

constant in variable a since there may be other vacua in fermion sectors. However there 
is a finite number of vacua hence G u{N \a,a 2 ) is at most a polynomial in a. We confirm 
this by explicitly calculating G u<yN \a,a 2 ) and G su<yN \a,a 2 ) for the lowest values of N. 
We have 
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G u ^(a,a 2 ) = l + a + a 2 + a\ G u ^ 3 \a,a 2 ) = 1 + a+a 2 + a 3 + a 5 + a 6 + a 7 -a 9 + a 10 + a n , 
and 

G su W(a, a 2 ) = l + a 2 -a 3 + a 4 , G su ^(a, a 2 ) = 1 + a 2 + a 5 + a 7 - a 8 + a 10 . 

The coefficients of polynomials G H (a, a 2 ) where H = U(N), SU(N) give us the difference 
between the number of bosonic and fermionic vacua. In general, since ,a 2 ) is a 

polynomial, the constraint for D^ B nF 's, coming from Eqn. ([1]), is now 

E (-i) B '^£,,», = o. 

for k greater then some k$. It also seems that 

E(-l) nF D H =01-1 

ng+2np=k 

for k < ko. 

We find it very interesting that although 's are just some group theory numbers, 

they are constrained by the dynamics of the properly chosen supersymmetric hamiltonian. 

3 The character method for D = 3 + 1,5 + 1,9 + 1 
SYMQM 

The generalization of the D = 2 case to D = 4, 6, 10 cases is now straightforward. The 
bosonic and fermionic creation operators a^, f^ a are now labeled by color index 6, spatial 
index i = 1, . . . , d and spinor index 

a = 1,2 forD=3+l, a = 1,2,3,4 for D=5+l , a = l,...,8 for D=9+l . 

The state with ub bosons and uf fermions is now 

n t bl n t 6 "B ft Cl ft C "F | M\ 

" H - • • " in B J oci ■ ■ ■ J a np I u /> 

and the number of U(N) singlets is 
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D U n i N lt = J d^ { K)XSy m {R n B B )XMt{R n F F )i 

where 

x(Rb) = d x d=1 (R), x(Rf) = (d-i) x d=l (R). 

The above equations for x(Rb) and x(Rf) ar e in fact the only difference between the 
d = 1 case. We can introduce the generating functions analogous to ([T]) and perform the 
same manipulations to find that the corresponding generating functions are 



From (|13p we identify the particle-hole symmetry 



Zi \d-l 



G u ^(a b) - (l^g^ r TT ^ TT(i - Zl ) { ^ (13) 



r,U(N),d _ n U(N),d 
^n B ,n F n B ,Ar 2 (d-l)-n F ' 



and supersymmetry 



V r>U(.N),d _ D U(N),d 

np even np odd 

Taking b = a and using the results from previous section we obtain 

G v W> d (a,a 



m v =i(i-« fc ) 

therefore 

np+ns=k 

The result does not depend on d which is surprising but possible since , a) is 

the generating function for the differences between bosonic and fermionic gauge invariant 
states. However, G u ( N ^' d (a,a) cannot be interpreted as the contribution to the Witten 
index for supersymmetric harmonic oscillator in d + 1 dimensions. This is because the 
supersymmetric harmonic oscillator with gauge degrees of freedom does not exist in d+1 > 
2 dimensions together with the singlet constraint. Such system cannot exist since the 
number of bosonic and fermionic degrees of freedom do not mach. For example there are 
d states with ns = 1 and np = 0, they are (a|), i = 1, . . . , d. On the other hand there 
are d — 1 states with = and np = 1, they are (fa), ol = 1, . . . , d — 1. The difference 
is precisely equal pjv(l) = 1- 
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riF 
n B 





1 


2 


3 


4 





1 


72 


28 


120 


406 


1 





288 


324 


2016 


9072 


2 


45 


3240 


3816 


21024 


89838 


3 


84 


12960 


23652 


150360 


692874 


4 


1035 


74520 


144000 


882720 


4049640 


5 


2772 


270864 


662436 


4331880 


20528802 


6 


16215 


1119096 


2906448 


18805104 


89459160 


7 


46530 


3635280 


10912572 


72993096 


353298330 


8 


189288 


12260160 


38914524 


259803720 


1263689658 



Table 1: Number of SU (2) singlets for D = 10 spacetime dimensions in sectors with < np < 4 
and < riB < 8. 

The generating function can be computed explicitly for arbitrary value of N although 
the general N dependence is difficult to obtain even for 6 = 0. The case N = 2 \s 
particularly easy to evaluate, we have 



dz 



1 d-2 



1 



d" 



2iri 



(zi - a) d (l - azx) 



2(d- 1)! dz° 



z d -\z-\f 
(1 - az) d 



(14) 



The cases with N = 3, 4, d = 3, 5, 9 are presented in Appendix E while the values of 



D 



SU (2),9 
ng,np 



for np < 12, riB < 10 are presented in Tables 1, 2 and 3. 



Some of these values have been already obtained earlier [T7] with considerable numer- 
ical effort by constructing the singlets directly . We see that the numbers of singlets grow 
extremely fast in this case, e.g. Dg 'j^' 9 ra 2.5 • 10 11 . It follows that the direct numerical 
approaches ( e.g. the cutoff method [7] ) to SYMQM in d = 9 dimensions is difficult to 
deal with even for the fastest computers. The d = 9 model is particularly troublesome be- 
cause the fermion number is not conserved hence one cannot diagonalize the hamiltonian 
in each fermion sector separately. 

If the hamiltonian has additional SO(d) symmetry then it is convenient to work with 
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Up 

n B 


5 


6 


7 


8 





1512 


4060 


8856 


17605 


1 


29232 


81648 


192528 


374544 


2 


321048 


907452 


2121192 


4230801 


3 


2426928 


6998040 


16742544 


33436080 


4 


14752080 


42942060 


103041000 


208064475 


5 


74701368 


220014900 


533991024 


1081967760 


6 


331885680 


984408096 


2399008272 


4891876599 


7 


1314510120 


3928885884 


9640642968 


19721394891 


8 


4754606472 


14285876220 


35181176976 


35181176976 



Table 2: Number of SU{2) singlets for D = 10 spacetime dimensions in sectors with 5 < np < 8 
and < n B < 8. 



Up 


9 


10 


11 


12 


E 


n B 















29512 


41392 


51520 


56056 


211068 


1 


626040 


908460 


1126944 


1205568 


4556448 


2 


7158600 


10328580 


12886776 


13896792 


51966252 


3 


56800008 


82741428 


103339320 


111140484 


414495042 


4 


355678200 


518416380 


649288080 


700074900 


2597347530 


5 


1856261448 


2718108792 


3408546960 


3673243476 


13592436138 


6 


8426386704 


12356524344 


15522106992 


16746545508 


61770199986 


7 


34066733976 


50095041876 


62995900968 


67964640282 


250166120589 


8 


72270195525 


125190973512 


184350316788 


232102914120 


920270148237 



Table 3: Number of SU{2) singlets for D = 10 spacetime dimensions in sectors with 9 < np < 
12 and < n B < 8. £ gives the cumulative size up to n B . 
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gauge and SO(d) singlets. The method presented here can be applied also in this case 
by computing the generating function for SU(N) invariant states with given spin ( see 
Appendix F ). Fortunately all the SYMQM have the SO{d) symmetry. It is then possible 
that the numerical analysis of these models is within our reach once we work in sectors 
with given SO (9) angular momentum. 



4 The partition functions for free Yang-Mills the- 
ories 

The generating functions G u ^' d {a, b) have an interesting application in, seemingly unre- 
lated, problem of computing the partition function of free Yang-Mills theories on S 1 xtime. 
We show in this section that the partition functions of such theories can be expressed in 
terms of G u ^' d (a, 6)'s in a rather simple way. Following Sundborg [20J and Aharony et 
al. [22] we write the partition function of the free Yang-Mills theory with n$, ny and n F 
number of scalar vector and fermion fields respectively as 

Z(x) = J d» G expi^^[Z B (x m ) + {-l) m Z F {x m )]x{R m )\ , (15) 

where Zb(x) and Z F (x) are bosonic and fermionic single partition functions given explic- 
itly, for D = 2, by 

Z B (x) = n s z s {x) + n v zv{x), Z F (x) = n F z F (x), 

1 + x 2 A^/x 

zs{x) = - , z v (x) = x , z F (x) = - • 

1 — X 1 — X 

It is clear that Z[x) is very similar to G u ^ N ^ ,d (a,b). The only difference is that there 

are no single particle partition function zb(x) and z F (x) in G u ^ N ^' d (a,b), instead there 

are two generating parameters a and b. In two dimensions scalar fields have the scaling 

dimension zero hance zs(0) = 1 and the partition function is divergent. To avoid this 

we take n$ = 0. Next, since z F (x) is a rational function, we Taylor expand z F (x) and 

substitute it to (j 151) . Using the formulas and conventions from Appendix A, (115H becomes 

N dz- z- °° (1 + x k+ ^ — ^ 4nr 



Z(x) 



li- ii mi 

- 1 j=l J i^j J k=0 i,3 K z o ' 



2 ■ 
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The above formula is more complicated then (|13j) due to the infinite product over k. 
However, the product does not appear when we take np = 0. In this case Z(x) becomes 



Z(x) =G u{N) > nv (x 2 ,0). (16) 

It is amusing that the partition function of free Yang-Mills theory on S 1 x R with ny 
vector fields is given directly by G u ^ ,nv (x 2 , 0). On the other hand two dimensional gauge 
theories have many exceptional features hence their partition functions may resemble some 
simplicity. Indeed, Yang-Mills theories on compact, orientable surfaces are exactly solvable 
[28j and their partition functions are known to be simple expressions depending on group 
theory parameters. 



5 Summary 

In this paper we focused on calculating the number of U(N) (SU(N)) singlets motivated 
by the numerical approach based on the cut Fock space. Once the basis is known the 
cutoff method may be used. However, the very construction of such a basis is far from easy 
and proved El [T7j to be very time consuming when symbolic programs are involved. 
The results presented here give us an algorithm to construct such basis theoretically 
thereby facilitating the numerical considerations. In particular we hope that the results 
presented here applied to the j = sector of D = 9 + 1 SYMQM will help to obtain the 
nonperturbative spectra of this highly complicated system. 

It is interesting that the generating functions G u ^ N ^' d (a,b) have other, seemingly in- 
dependent, applications, i.e. they give rise to the Witten index in a class of models as 
well as they can be used to compute the partition functions of free Yang- Mills theories on 
S 1 x R. 
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7 Appendix A - Group theory conventions and 
the integral representation of G u ( N \a,b) 

Here we give the conventions used in integrals over characters and derive ([2]). They can 

be found in, e.g. ref. [29]. The methods used in this appendix are similar to the ones 

used in [22] . 

The U(N) invariant measure is 

1 N da I N da N 

dVu(N) = J^ U^ I M I'' d »SU(N) = V^E"*) I M | 2 , a t e [0,2vr], 

i=l i=l i=l 

where 5p is a periodic Dirac delta with period 2tt 

S P (x) = ^<5(a;-27rA;), 

k 

the measure factor M is given by Vandermonde determinant 



M = Det(zf i} ) = Y[(zi - Zj ), Zj 

i<j 



the symmetric and antisymmetric powers of R, Xslm(^) an d xJ^C^)) are given by 
Frobenius formulas 

3%. v TT 1 xlk{Rk) 

X S ym{ R )- X. 11 "I k i k ' 
J2 k ki k =n B k=l 

and the characters \ are given by Weyl determinant formula 

Det(z {N - i+h) ) 
X (R) = X (Wi}lt) = ! (N .J , X ^ = 

The numbers enumerate the representation in which the character is calculated. In 
our case it is the adjoint representation of U(N) ( or SU(N) ) therefore h, ■ ■ ■ , In) = 
(2,1,..., 1, 0) . In this representation the characters simplify into 



XU(N)({Ui}) = E~' 



*,3 
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Xsu(N)({ai}) = ~ L 

. . Zj 

hi 

In order to derive ([2|) let us introduce 

oo 

F Sym {a,{ai}tx) = anB XSym(R nB ), \a\<l, 



n B =0 



np=0 

G H (a,b)= / dfi G F S ym(a,{ai})F A it(b,{ai}), 

J[0,2n] N 

D n = 1 (-1)"* ^ Gf 

The last sum is in fact finite since for U(N) ( or SU(N) ), XAit(R nF ) = when 
rip > N 2 ( or nj7 > iV 2 — 1). It is however more convenient to work with infinite sum as 
we will see in the following. The b variable is not bounded. 

Using the standard manipulations we obtain 

/ 00 a k \ 
F S ym(a,{ai}f =1 ) = exp |^ — X {{kai}f =1 ) J , 

F M t{b,{a.i}f =l ) = — N ■ 

The generating function can be calculated explicitly for arbitrary U(N) and SU(N). We 
have 



F, 



nu(i-«f 



2i \ ' 



Fsyl N \a, = (1 - a)*ffi?(o, {a i} f =l)j 



therefore the generating functions are 



7 



i^j ■> K z j' 
where in the last integral we changed variables z% — ► z%j rX?=i z j-> Z N — ► 11^=1 z j- 
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8 Appendix B - Explicit calculation of G u ( N \a,0) 



Here we evaluate the integral ^fy explicitly for 6 = 0. It can be done with use of the 



Cauchy determinant formula 
del 

which for U (N) gives 



N(N-l) „ 

1 ) = a 2 TT-TT {Zl ~ Zj) , (17) 

Zi - azj J (a - l) N 11 Zi 11 (zi - azj)(zj - aZi) ' 



G UW (a,0) = -Li=£> / n^tot—L-). (18) 



a 2 

The determinant under the integral (j!8f) can be expressed as a sum over cycles. The 
integration over each cycle can be done separately and it gives the factor fc , i.e. 



J \ZS =1 .-1 



dzj 1 



*.i=i M 27ri z\ - az 2 z 2 - az 3 " ' Zk-i - azk Zk - az\ 1 - a k ' 
therefore we obtain 



where Li lt „i N is the number of different permutations with the same cycle structure given 
by the partition (T 1 . . .N iN ) i.e. L^„ An = iVI/tT^Li The right hand side of {TSJ 

is in fact very simple due to the Cayley identity £ 



N i N 1 

(-1) E — II 773^ = (-l^a-^JV! II TTT^y (20) 

ji+2i 2 + ...+A r «„=A r fe=l V 7 fe=l V ; 

It can be proven most efficiently with use of the Bell polynomials |30j . Therefore we finally 
obtain @. 



5 The general form of Cauchy determinant formula is 

, , , 1 \ Ui<j( z i ~ z j)i x i ~ x o) 
det ' ' — 



Zi Xj J Wi^(%i x j) 

which for X4 = azi yields (|17[) . 

6 The Cayley identity by definition is 

N 1 N I 

e £ <i-^ n (i_ a fe )4fc =iv! n (i_ a fc)> 

ii+2i 2 +---+iVi„=A f fc=l V ' fe=l V 7 

which for a — ► 1/a yields ([20]) . 
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For SU(N) the only difference is that Tr(A) = hence in Eqn.(|3|) there is no 1/(1 — a) 
factor, i.e. 

oo 

G su W(a,0) = (1 - a)G u ^(a,0) = £ a n *q N (n B ), 

n B =0 

where ^(n^) is the number of partitions of ns into numbers 2, 3, . . . , N. 
The generating function G su ^ (a, b) clearly have the form 

/ N \ N 2 -l 

\fc=2 V V i=0 

where c^ U ^ N \a) are polynomials in variable a. 



9 Appendix C - Examples of polynomials c^ 7 ^ 

and c su ^ 

Here we list the polynomials and c su ^ for = 2,3,4. They can be obtained 

from equation ([2]) using some symbolic program, e.g. Mathematica, to evaluate the cor- 
responding residues. For N = 2 they are 

c£ (2) = l, c^ = l + a, c U 2 W=2a, cp> = c£f '. 



and 



For Af = 3 they are 



SU{2) _ 1 5C/(2) _ 5C7(2) _ St/(2) 



c ^ 3 ) = 1, % P) = 1 + a + a 2 j ^(3) = 2fl + 2a 2 + 2a3) 

= 1 + 2a + 3a 2 + 5a 3 + a 4 , c 4 ^ 3) = 1 + 3a + 6a 2 + 5a 3 + 3a 4 , cf (3) = cfi? . 



and 

c 5E/(3) = ^ ^(3) = Q + ^(3) = a + a 2 + 2a3) 

c f (3) = ! + a + 2a 2 + 3a 3 + a 4 ; ^(3) =2a + 4fl 2 + 2fl 3 + 2fl 4 ^(3) = ^(3) 

For A^ = 4 they are 

c ^ 4 ) = l, c f W = l + a + a 2 + a 3 , c^ (4) = 2a + 2a 2 + 4a 3 + 2a 3 + 2a 5 , 
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^ (4) = i + 2a+4a 2 +8a 3 +8a 4 +7a 5 +5a 6 +a 7 , c^ (4) = l+3a+9a 2 +13a 3 +19a 4 +17a 5 +18a 6 +7a 7 +3a 8 , 



1 + 4a + 11a 2 + 22a 3 + 33a 4 + 38a 5 + 34a 6 + 23a 7 + 11a 8 + 3a 9 



c^ (4) = 1 + 5a + 12a 2 + 33a 3 + 45a 4 + 62a 5 + 55a 6 + 45a 7 + 22a 8 + 11a 9 + a 10 , 
c^ {4) = 1 + 5a + 16a 2 + 37a 3 + 59a 4 + 75a 5 + 77a 6 + 60a 7 + 37a 8 + 17a 9 + 4a 10 , 
c^ {4) = 2 + 4a + 18a 2 + 36a 3 + 68a 4 + 78a 5 + 86a 6 + 64a 7 + 46a 8 + 18a 9 + 6a 10 . 

and 

cf ( 4 ) = 1, c f W = a + a 2 + a 3 , cf (4) = a + a 2 + 3a 3 + 2a 4 + 2a 5 , 
cf (4) = i +0+ 3 a 2 + 5 a 3 + 6a 4 +5a 5 +5a 6 +a 7 , cf (4) = 2a+6a 2 +8a 3 +13a 4 +12a 5 +13a 6 +6a 7 +3a s 

5J7(4) = 1 + 2a + 5a 2 + Ua 3 + 2Qa 4 + 2 g a 5 + 21fl 6 + 17ffl 7 + 8a 8 + g fl 9 



3a + 7a 2 + 19a 3 + 25a 4 + 36a 5 + 34a 6 + 28a 7 + 14a 8 + 8a 9 + a 10 , 



c SU(4) _ 1 + 2a + g a 2 + 18a 3 + 34a 4 + 39a 5 + 43fl 6 + 32q 7 + 23a 8 + g a 9 + 3a 10_ 



10 Appendix D - Linear dependence of the trace 
operators for the SU(3) group 

Here we derive the linear dependence of 17 trace operators listed in subsection 2.1.1. Our 
starting point is the equation 



aV/T + /W + at/V = ^a f )f + (aW + (aV/t). (21) 
It is the source of the following relations. 



Multiplying (|2ip from the right hand side by f* f> and taking the trace gives 
(at/tat/t/t) = I( a tat)(/t/t/t) + ( t/t)( t/t/t) - 2(a+a+/+/t/t). 
Therefore, we may neglect, e.g. (a^/^a^/^/^)- 



Multiplying (|21f) from the right hand side by /V T an d from the left hand side by a) and 
then taking the trace gives 

(atat/tat/t/t) + (at/tatat/t/t) = (at/t)(atat/t/t) + ( a tat/t)(at/t/t)-I( a t a t a t )(/ t / t / t )i 
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where the Cayley-Hamilton theorem for a' matrices was also used. 

Multiplying (I21D from the right hand side by f'a'f' and taking the trace gives 
(«t a t / t / t a t / t )+(a t / t a t / t a t / t )+(a t a t / t a t / t / t ) = I( a t a t )(a t / t / t / t ) _ (a t a t / t )(a t / t / t ) . 

From the two above equations it follows that (at/tat /t a t/t) can be expressed in terms of 
multiple trace operators. Also, we may neglect, e.g. (a'a* f'w f' f'). 

Multiplying (|2T1) from the right hand side by a) a) f> f> and taking the trace gives 
(at a t/t a tat/t/t) = i(at/t)(at a t)(at/t/t) + (At/t^aW) - |(ata+a+)( a t/t/t/t) 

-I( a ta+)(at/t a t/t/t), 

where the Cayley-Hamilton theorem for at matrices was also used. Therefore we neglect 
(atat/tatat/t/t). 

Multiplying (|2ip from the right hand side by /t a t a t/t a nd taking the trace gives 
2(atat/tat fl t/t/t) + (at a t/taW/t) = ±(<J<J)(<J<J f^) + i(at/t)( a t/t a +at/t) 

-(atat/t)(atat/t/t), 

where the Cayley-Hamilton theorem for at matrices was again used. From the two above 
equations it follows that we can also neglect (at at/tat /tat/t). 

Multiplying (|21f) from the right hand side by at /tat/t and taking the trace gives 
2(ataW/t t/t) + i(atat)(/t a t/t a t/t) + i(atat t)(/t/t a t/t) = i( a t a t)(/t a t/tat/t) 

+(at/t)(atat/t a t/t), 

where the Cayley-Hamilton theorem for at matrices was used. The operator (at at/tat /tat/t) 
was already excluded before hence, from the above identity, it follows that there is 
a relation between the multiple trace operators. Therefore we can neglect one, e.g. 
(ot/t)(otot/tot/t). 

Multiplying (|2ip from the right hand side by a/cft /tat/t and taking the trace gives 
(atat/tatat/tat/t) = -I(at/t)(atatat)(at/t/t)+(atat/t)(atat/tat/t)-^(atatat)(at/tat/t/t) 
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--(atat)(aW/W), 

where the Cayley-Hamilton theorem for a) matrices was used. Therefore, we neglect the 
operator (aW/taW/W/t). 

Multiplying fl21f) from the right hand side by a) f'a'a* /' and taking the trace gives 
2(atat/ta+at/tat/t) + i(atat)(/tat/tatat/t)+i(atatat)(/t/t a tat/t) = i( a t a t )(/ t a t / t a t a t / t ) 

+(atat/t)(ot/t +at/t), 

where the Cayley-Hamilton theorem for a) matrices was used. The (a' a) a) p p) is 
already excluded therefore we can neglect one multitrace operator, e.g. (atoV)(^<i'/ s / )■ 

Multiplying (f2Tj) from the right hand side by a)a) pa/w f> and taking the trace gives 
(aV/W/W/t) + ^(a t ot)(aV t a t / t a t a t / t ) + ~(aW)(at/t/W/t) 

+i(otot)(/toto+/taV/t) + i(atatat)(/t a t/totot/+) = ±(<J<J)(fWfWfi) 
+ ±( a t/t)( a t a t)( a t / t a t a t / t) + i( a t/t)( a t a t a t) (/ t a t a t / t )) 

where the Cayley-Hamilton theorem for at matrices was also used. It follows that the 
operator (a^a^f^a^a'pa'a'f') can be excluded as well. 

There are no other independent relations following from (|21|) except from the ones 
where the Cayley-Hamilton theorem for a) matrices is used. We excluded 9 operators as 
being linearly dependent. The remaining 8 trace operators are indicated in subsection 
2.1.1 . 



11 Appendix E - Examples of generating func- 
tions G SXJ ^ d (a,b) 

Below we list the generating functions G su ( N ^' d (a,b), for d = 3, 5,9, N = 2,3. The case 
of U(N) group is obtained from the identity 

G u W> d (a,b) = ^f^- G su ^' d (a,b). 
(1 a) 
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For d = 3, N = 2 we have 

G SU(2),d=3( b) = 1 S^ b i SU(2),d=3 

, SU(2),d=3 SU(2),d=3 l 

wnere = c 6 _> and 

^(2)^=3 = 1 _ a + fl2) ^(2)^=3 = c 5C/(2),,=3 = : + 8fl + 7fl 2 + 2fl 3 _ 2fl 4 _ ^ 

c 5C/(2),d=3 = 4 + 2a + 16a 2 + 4a 3 _ 4a 4 _ 2fl 5 

For d = 5, N = 2 we have 



We will not list the rest of Cj's since there are many of them and they become more 
complicated. 

For d = 9, N = 2 we have 

24 

u ^ b > (i _ a yn _ Q 2)i7 2^ ^ °> 

v ; v ; i=0 



SC/(2),d-9 = 1 _ 7a + 49a 2 _ M7a 3 + 441a 4 _ 735a 5 + 1225fl 6 _ 1225fl 7 + +12 25a 8 - 735a 9 



+441a 10 - 147a 11 + 49a 12 - 7a 13 + a 14 . 
For d = 3, TV = 3 we have 

G SU(3),d=3 ( b) = 1 Vf-lV c SU (^ d = s b i 

c 5C/(3),d=3 = 1 _ a _ 2a 2 + ga 3 + g a 4 _ g a 5 + fl 6 + 1?a 7 + fl 8 _ g a 9 + Q a W + 6a ll _ 2ffl 12 _ fl 13 + a 14 _ 

For d = 5, N = 3 we have 

32 

U ^ b > (l_ a )3 (1 _ a 2)16 (1 _ a 3)13Z^^ L ) C i °> 

c SU(3),d=5 = 1 _ 3a+ 2 a 2 + 3 4a 3_ 4a 4_ 18a 5 +421a 6 +624a 7 +251a 8 +2107a 9 +5377a 10 +4766a ll 

+6384a 12 +16031a 13 +19327a 14 +14592a 15 +21381a 16 +29839a 17 +21381a 18 +14592a 19 +19327a 20 
+ 16031a 21 +6384a 22 +4766a 23 +5377a 24 +2 107a 25 +251a 26 +624a 27 +421a 28 -18a 29 -4a 30 +34a 31 
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+2a 32_3a 33 + a 34 . 

For d = 9, N = 3 we have 

1 64 

G SU(3),d=9 (a b) = i ^ ( y SU(3),d=9 b i 

where c ^ u ^' d ~ 9 [ s Q f order 74. 

12 Appendix F - Number of gauge singlets with 
given angular momentum 

Here we discuss the character method applied to sectors with fixed angular momentum. 
The projection to sectors with fixed angular momentum j is due to the decomposition 

V = Symfr^Af 1 ) x Alt{®£i*! =1/2 ), 

where Ajr 1 , Fj'~ 1 ^ 2 are vector spaces spanned by ajr 1 | 0), and // _1 ^ 2 | 0) where operators 
°fc _1 ' // _1//2 are assumed to carry 50(3) spin 1 and 1/2 respectively. Therefore the 
dimensions of subspaces with angular momentum j are 

where dnsu(N) an d dfisotd) are SU(N) and SO(d) invariant measures. We will restrict 
to the d = 3 case hence we take 

1/3 f 
d Vso(3) = ~ sin2 g d & ^ G t ' 27r ]' / ^50(3) = !■ 

i?-^ -1 , R^ 1 ^ 2 are the adjoint representation of SU(N) and j = 1, j = 1/2 representations 
of SO(d) respectively, i.e. 

x(^ =1 ) = x(R so{d) > j=1 )x(R sum J=1 ), 
x(i?]T 1/2 ) = x(R so(d) > ^ 2 )x(R sum J=1 ). 

For d = 3 we have 

sm^a / — 
2 fc=-j 



27 



therefore for SU (2) gauge group we obtain 

X (4 =1 ) = x(R SO(3) ' J=1 )X(R SU{2) > j=1 ) = (l + 2cos/?)(l + 2cosa), 

X (4 =1 / 2 ) = X (R so ^> i=V*) x ( R Sum, i=i) = 2cos^(l + 2cosa). 

The explicit calculation of the generating function for Dfj„,J rf_3, J is now straight- 
forward. Below we perform the calculation for the purely bosonic sector, i.e. we evaluate 

oo oo 

Tlfl=0 j = 

Using the same conventions and techniques as in Appendix A we perform the sum over 
rig. We have 

k=j 



G(a,c) =XXy d f 1 so(3)d^su(2) ( I F ( a ' z; c ' *)> 



z = e ia , t = e 



l3 



3=0 \k=-j 

where 

F(a, z; c, t) 



t 3 z 3 



(1 — a)(l — at)(l — az)(l — atz)(a — t)(a — z)(a — tz)(at — z)(t — az) 
Now the sum over j is also possible and the evaluation of the resulting integral gives 



G(a,c) 



1 - a 2 c + a 4 c 2 



(1 - a 2 )(l - a 3 )(l - a 4 )(l - a 2 c)(l - a 2 c 2 ) ' 
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